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Abstract With the recent availability of state-of-the-art heavy-ion stable and ra- 
dioactive beams, there has been a renew interest in the investigation of nuclear re- 
actions with heavy ions. I first present the role of inelastic and transfer channel cou- 
plings in fusion reactions induced by stable heavy ions. Analysis of experimental 
fusion cross sections by using standard coupled-channel calculations is discussed. 
The role of multi-neutron transfer is investigated in the fusion process below the 
Coulomb barrier by analyzing 32 S+ 90 96 Zr as benchmark reactions. The enhance- 
ment of fusion cross sections for 32 S+ 96 Zr is well reproduced at sub-barrier energies 
by NTFus code calculations including the coupling of the neutron-transfer channels 
following the Zagrebaev semi-classical model. Similar effects for 40 Ca+ 90 Zr and 
40 Ca+ 96 Zr fusion excitation functions are found. The breakup coupling in both the 
elastic scattering and in the fusion process induced by weakly bound stable projec- 
tiles is also shown to be crucial. In this lecture, full coupled-channel calculations 
of the fusion excitation functions are performed by using the breakup coupling for 
the more neutron-rich reaction and for the more weakly bound projectiles. I clearly 
demonstrate that Continuum-Discretized Coupled-Channel calculations are capable 
to reproduce the fusion enhancement from the breakup coupling in 6 Li+ 59 Co. 



1 Introduction 

Heavy-ion fusion reactions at bombarding energies at the vicinity and below the 
Coulomb barrier have been widely studied JT] |2] [3] @] |5] . In low-energy fusion reac- 
tions, the very simple one-dimensional barrier-penetration model (1D-BPM) IHEl 
is based upon a real potential barrier resulting from the attractive nuclear and repul- 
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sive Coulomb interactions. For light- and medium-mass nuclei, one only assumes 
that the di-nuclear system (DNS) fuses as soon as it has reached the region inside the 
barrier i.e. within the potential pocket. If the system can evolve with a bombarding 
energy high enough to pass through the barrier and to reach this pocket with a rea- 
sonable amount of energy, the fusion process will occur after a complete amalgation 
of the colliding nuclei forming the compound nucleus (CN). On the other hand, for 
sub-barrier energies the DNS has not enough energy to pass through the barrier. 

In reactions induced by stable beams, the specific role of multi-step nucleon- 
transfers in sub-barrier fusion enhancement still needs to be investigated in detail 
both experimentally and theoretically fl6] |7] [8] |9] [10] [TT] [121 [T3l [14]| . In a complete 
description of the fusion dynamics the transfer channels in standard coupled-channel 
(CC) calculations J2] [8j [TO] [T4] \15l have to be taken into account accurately. It is 
known, for instance, that neutron transfers may induce a neck region of nuclear 
matter in-between the interacting nuclei favoring the fusion process to occur. In this 
case, neutron pick-up processes can occur when the nuclei are close enough to inter- 
act each other significantly ]7][8), if the Q-values of neutron transfers are positive. It 
was shown that sequential neutron transfers can lead to the broad distributions char- 
acteristic of many experimental fusion cross sections. Finite Q-value effects can lead 
to neutron flow and a build up of a neck between the target and projectile JS). The 
situation of this neck formation of neutron matter between the two colliding nuclei 
could be considered as a "doorway state" to fusion. In a basic view, this intermediate 
state induced a barrier lowering. As a consequence, it will favor the fusion process 
at sub-barrier energies and enhance significantly the fusion cross sections. Experi- 
mental results have already shown such enhancement of the sub-barrier fusion cross 
sections due to the neutron-transfer channels with positive Q-values J6]|9l. 

In reactions induced by weakly bound nuclei and/or by halo nuclei, the influ- 
ence on the fusion process of coupling both to collective degrees of freedom and to 
transfer/breakup channels is a key point J5] H] [5] for the understanding of N-body 
systems in quantum dynamics 1 1 1 . Due to their very weak binding energies, a diffuse 
cloud of neutrons for 6 He or an extended spatial distribution for the loosely bound 
proton in 8 B would lead to larger total reaction (and fusion) cross sections at sub- 
barrier energies as compared to 1D-BPM model predictions. This enhancement is 
well understood in terms of the dynamical processes arising from strong couplings 
to collective inelastic excitations of the target (such as "normal" quadrupole and oc- 
tupole modes) and projectile (such as soft dipole resonances). However, in the case 
of reactions where at least one of the colliding nuclei has a sufficiently low binding 
energy for breakup to become a competitive process, conflicting conclusions were 
reported UUffllU El H3. 

Recent studies with Radioactive Ion Beams (RIB) indicate that the halo nature 
of 6 8 He |[T8l[T9ll20ll2Tll22ll . for instance, does not enhance the fusion probability 
as anticipated. Rather the prominent role of one- and two-neutron transfers in 6,8 He 
induced fusion reactions was definitively demonstrated. On the other hand, the effect 
of non-conventional transfer/stripping processes appears to be less significant for 
stable weakly bound projectiles. Several experiments involving 9 Be, 7 Li, and 6 Li 
projectiles on medium-mass targets have been undertaken. 
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Fig. 1 Comparison be- 
tween the fusion-evaporation 
(ER) excitation functions of 
32 S+ 90 Zr (open circles) and 
32 S+ % Zr (points) as a func- 
tion of the center-of-mass 
energy. The error barrs of the 
experimental data taken from 
Ref. 1231 represent purely 
statistics uncertainties. (Cour- 
tesy of H.Q. Zhang) 
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Fig. 2 Ratios of measured 
fusion cross sections for 6 Li 
and 7 Li projectiles with 24 Mg, 
28 Si and 59 Co targets as a 
function of E c m ./Vj,. The 
solid line gives the 1D-BPM 
prediction while the dotted 
line shows results obtained 
from Wong's prescription. 
(This figure originally shown 
in Ref. (24) for 6 7 Li+ 59 Co has 
been adapted to display com- 
parisons with other lighter 
targets |25l|26l|27l|28l) 
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2 Experimental results 

In this lecture we first present the role of inelastic and transfer channel couplings 
in experimental data obtained in fusion reactions induced by stable 32 S projectiles 
l23l . The breakup coupling in both elastic scattering data and in the fusion data are 
also shown for weakly bound 6 7 Li projectiles l24l . 



2.1 i2 S + w Zr and i2 S + 9b Zr reactions 

In order to investigate the role of neutron transfers we further study 32 S + 90 Zr and 
32 S + 96 Zr as benchmark reactions. Fig. 1 displays the measured fusion cross sec- 
tions for 32 S + 90 Zr (open circles) and 32 S + 96 Zr (points). We present the analysis of 
excitation functions of evaporation residues (ER) cross sections recently measured 
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with high precision (i.e. with small energy steps and good statistical accuracy for 
these reactions [23]). 

The differential cross sections of quasi-elastic scattering (QEL) at backward an- 
gles were previously measured by the CIAE group ifPTl . The analysis of the corre- 
sponding BD-QEL barrier distributions (see solid points in Fig. 3) already indicated 
the significant role played by neutron tranfers in the fusion processes. 

In Fig. 3 we introduce the experimental fusion-barrier (BD-Fusion) distributions 
(see open poins) obtained for the two reactions by using the three-point difference 
method of Ref. [8] as applied to the data points of Ref. 12311 plotted in Fig. 1 . It is in- 
teresting to note that in both cases the BD-Fusion and BD-QEL barrier distributions 
are almost identical up to E c . m . ~ 85 MeV. 



2.2 6 Li + 59 Co and 7 Li + 59 Co reactions 

The fusion excitation functions were measured for the 6 ' 7 Li+ 59 Co reactions 11241 at 
the VIVITRON facility of the IPHC Strasbourg and the Pelletron facility of Sa5 
Paulo by using y-ray techniques. Their ratios are presented in Fig. 2 with compar- 
isons with other lighter targets |25] |26] |27] |28] . The theoretical curves (1D-BPM 
IH 12 and Wong ifTSI do not take into account the breakup channel coupling that is 
discussed in one of the following sections in more details. 



3 Coupled channel analysis 

Analysis of experimental fusion cross sections by using standard CC calculations is 
first discussed with the emphasis of the role of multi-neutron transfer in the fusion 
process below the Coulomb barrier for 32 S+ 90 96 Zr as benchmark reactions. 



3.1 i2 S + 90 Zr and i2 S + 9b Zr reactions 

A new CC computer code named NTFus |29l taking the neutron transfer channels 
into account in the framework of the semiclassical model of Zagrebaev iflOl has been 
developed. The effect of the neutron transfer channels yields a fairly good agreement 
with the data of sub-barrier fusion cross sections measured for 32 S + 96 Zr, the more 
neutron-rich reaction 11231 . This was initially expected from the positive Q-values 
of the neutron transfers as well as from the failure of standard CC calculation of 
quasi-elastic barrier distributions without neutron-transfers coupling lf]~3l as shown 
by the solid line in Fig. 3(b). 

By fitting the experimental fusion excitation function displayed in Fig. 1 with 
NTFus CC calculation ll29l . we concluded l30l that the effect of the neutron transfer 
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channels produces significant enhancement of the sub-barrier fusion cross sections 
of 32 S + 96 Zr as compared to 32 S + 90 Zr. A detailed inspection of the 32 S + 90 Zr fu- 
sion data presented in Fig. 1 along with the negative Q-values of their corresponding 
neutron transfer channels lead us to speculate with the absence of a neutron trans- 
fer effect on the sub-barrier fusion for this reaction. With the semiclassical model 
developed by Zagrebaev ifTDl we propose to definitively demonstrate the significant 
role of neutron transfers for the 32 S + 96 Zr fusion reaction by fitting its experimental 
excitation function with NTFUS code ll29l calculations, as shown in Fig. 3. 



Fig. 3 Barrier distributions 
(BD) from the fusion ER 
(open circles) cross sections 
1231 . plotted in Fig.l, and 
quasielastic scattering (solid 
circles) cross sections 1131 
for 32 S+ 90 Zr (a) and 32 S+ 96 Zr 
(b). The dashed and solid 
black lines represent uncou- 
pled calculations (1D-BPM) 
and the CC calculations with- 
out neutron transfer coupling. 
The red dash-dotted line rep- 
resents the CC calculations 
with neutron transfer coupling 
for the 32 S+ % Zr reaction. 
(Courtesy of H.Q. Zhang) 
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The new oriented object NTFUS code ll29l , using the Zagrebaev model ifTUl was 
implemented (at the CIAE) in C++, using the compiler of ROOT 11311 . following 
the basic equations of Ref. l32l . Let us first remind the values chosen for the de- 
formation parameters and the excitation energies that are given in Refs. EI [33l [34l 
(see Tables given in |30| for more details). The quadrupole vibrations of both the 
90 Zr and 96 Zr are weak in energy; they lie at comparable energies. The 96 Zr nucleus 
presents a complicated situation l35l : its low-energy spectrum is dominated by a 2 + 
state at 1.748 MeV and by a very collective [B(E3;3~ ->■ 0+) = 51 W.u.] 3~ state 
at 1.897 MeV. CC calculations explained the larger sub-barrier enhancement as due 
mainly to the strong octupole vibration of the 3~ state in 36 S + 96 Zr 11361 . However, 
the agreement is not so satisfactory below the barrier for 32 S + 96 Zr (see solid line 
of Fig. 3.b), as well as for 40 Ca + 96 Zr |9j and, therefore, there is the need to take 
neutron transfers into account. 

The main functions of the code NTFus are designed to calculate the fusion ex- 
citation functions with normalized barrier distribution (based on experimental data) 
given by CCFULL lfT5l , we take the dynamical deformations into account. In order 
to introduce the role of neutron transfers, the NTFus code |29l applies the Za- 
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grebaev model IflOl to calculate the fusion cross sections Of us {E) as a function of 
center-of-mass energy E. Then the fusion excitation function can be derived using 
the following formula ifTUl : 

T/(E) = 

r 1 _ rQo(k) 
fWjrL a k {E,l,Q)xP HW (B,E + Q,l)dQdB. (1) 

and 

O fus (E) = ?—Z(2l + m(E). (2) 

where T/(E) are the transmission coefficients, E is the energy given in the center- 
of-mass system, B and f(B) are the barrier height and the normalized barrier distri- 
bution function, Phw is the usual Hill- Wheeler formula. / is the angular momentum 
whereas l cr is the critical angular momentum as calculated by assuming no coupling 
(well above the barrier). Ok(E,l,Q) and Qo{k) are, respectively, the probabilities 
and the Q-values for the transfers of k neutrons. And 1 /N tr is the normalization of 
the total probability taking into account the neutron transfers. 

The NTFus code |29l uses the ion-ion potential between two deformed nuclei 
as developped by Zagrebaev and Samarin in Ref. [32|. Either the standard Woods- 
Saxon form of the nuclear potential or a proximity potential |37|) can be chosen. 
The code is also able to predict fusion cross sections for reactions induced by halo 
projectiles [30); for instance 6 He + 64 Zn l22l[38l . In the following, only comparisons 
for 32 S + 90 Zr and 32 S + 96 Zr are discussed. 

For the high-energy part of the 32 S + 90 Zr excitation function, one can notice 
a small over-estimation of the fusion cross sections at energies above the barrier 
up to the point used to calculate the critical angular momentum. This behavior can 
be observed at rather high incident energies - i.e. between about 82 MeV and 90 
MeV (shown as the dashed line in Fig. 3. (a) for 32 S + 90 Zr reaction). We want to 
stress that the corrections do not affect our conclusions that the transfer channels 
have a predominant role below the barrier for 32 S + 96 Zr reaction, as shown by the 
dotted-dashed red curve in Fig. 3.(b). 

As expected, we obtain a good agreement with calculations not taking any neu- 
tron transfer coupling into account for 32 S + 90 Zr as shown by the solid line of 
Fig. 3. (a) (the dashed line are the results of calculations performed without any cou- 
pling). On the other hand, there is no significant over-estimation at sub-barrier ener- 
gies. As a consequence, it is possible to observe the strong effect of neutron transfers 
on the fusion for the 32 S + 96 Zr reaction at sub-barrier energies. Moreover, the bar- 
rier distribution function f(B) extracted from the data contains the information of 
the neutron transfers. These information are also contained in the transmission co- 
efficients, which are the most important parameters for the fusion cross sections to 
be calculated accurately. The f(B) function as calculated with the three-point for- 
mula [8 1 will mimic the differences induced by the neutron transfer taking place in 
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sub-barrier energies where the cross section variations are very small (only visible if 
a logarithm scale is employed for the fusion excitation function). It is interesting to 
note that the Zagrebaev model HI 01 implies a modification of the Hill- Wheeler prob- 
ability and does not concern the barrier distribution function f(B). Finally, the code 
allows us to perform each calculation by taking the neutron transfers into account 
or not. 

The calculation with the neutron transfer effect is performed up to the channel 
+4n (k=4), but we have seen that we obtain the same overall agreement with data up 
to channels +5n and +6n (30l . As we can see on Fig. 3.(b), the solid line representing 
standard CC calculations without the neutron transfer coupling (the dotted line is 
given for uncoupled calculations) does not fit the experimental data well at sub- 
barrier energies. On the other hand, the dotted line displaying NTFus calculations 
taking the neutron transfer coupling into account agrees perfectly well with the data. 
As expected, the Zagrebaev semiclassical model's correction applied at sub-barrier 
energies enhances the calcutated cross sections. Moreover, it allows to fit the data 
reasonably well and therefore illustrates the strong effect of neutron transfers for the 
fusion of 32 S + 96 Zr at subbarrier energies. 

The present full CC analysis of 32 S + 96 Zr fusion data |23][30) using NTFus (29) 
confirms perfectly well first previous CC calculations lITOl describing well the earlier 
40 Ca + 90 ' 96 Zr fusion data [9] and, secondly, very recent fragment- 7 coincidences 
measured for 40 Ca + 96 Zr multi-neutron transfer channels 11351 . 



3.2 6 Li + 59 Co and n Li + 59 Co reactions 

For reactions induced by weakly bound nuclei lfT6l rT7l |25l |26l |27l |28l [39l |40l |4T1 
and exotic nuclei HU |2D1 |2U HU EU SI E3 ggl , the breakup channel is 
open and plays a key role in the fusion process near the Coulomb barrier similarly 
to the transfer-channel coupling described in the previous section. It is therefore 
appropriate to use the Continuum-Discretized Coupled-Channel (CDCC) approach 
[47 48 49 50 1 to describe the influence of the breakup channel in both the elastic 
scattering and the fusion process at sub-barrier energies. 

Theoretical calculations (including CDCC predictions given in Refs. l47l l49l 
indicate only a small enhancement of total fusion for the more weakly bound 6 Li 
below the Coulomb barrier (see curves of Fig. 2), with similar cross sections for both 
6 ' 7 Li+ 59 Co reactions at and above the barrier 11241 . It is interesting to notice, how- 
ever, that the same conclusions have been reached for other targets such as 24 Mg 
11261 and 28 Si l|25l |27l |28l as can be clearly seen in the plot of Fig. 2. Thess results 
are consistent with rather low breakup cross sections measured for the 6,7 Li+ 59 Co 
reactions even at incident energies larger than the Coulomb barrier [ 39, 40 41 1. But 
the coupling of the breakup channel is extremely important for the CDCC analy- 
sis of the angular distributions of the elastic scattering ll47l as shown in Fig. 4 for 
6 Li+ 59 Co. The curves show the results of calculations with (solid lines) and with- 
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Fig. 4 Ratios of the elastic scattering cross-sections to the Rutherford cross sections as a function 
of cm. ana le for the 6 Li+ 59 Co system (47]. The curves correspond to CDCC calculations with 
(solid lines) and without (dashed lines) 6 Li — > a + breakup couplings to the continuum for 
incident 6 Li energies of (a) 30 MeV, (b) 26 MeV, (c) 18 MeV and (d) 12 MeV. (This figure has 
been adapted from the work of Ref. [47 1) 



out (dashed lines) 6-7 Li — > a + d,t breakup couplings. The main conclusion is that 
effect of breakup on the elastic scattering is stronger for 6 Li than 7 Li. 

A more detailed investigation of the breakup process in the 6 Li+ 59 Co reaction 
with particle coincidence techniques is now proposed to discuss the interplay of fu- 
sion and breakup processes. Coincidence data compared to three-body kinematics 
calculations reveal a way how to disentangle the contributions of breakup, incom- 
plete fusion and/or transfer-reemission processes ll39] l40l |4T1 . 

Fig. 5 displays experimental (full rectangles) and theoretical angular distributions 
(solid lines) for the sequential (SBU) and direct (DBU) projectile breakup processes 
at the two indicated bombarding energies for the 6 Li+ 59 Co reaction. In the CDCC 
calculations the a + d binning scheme is appropriately altered to accord exactly with 
the measured continuum excitation energy ranges. For this reaction it was not neces- 
sary to use a sophisticated four-body CDCC framework. The CDCC cross sections 
l47l are in agreement with the experimental ones llT6l |40] |4T1 , both in shapes and 
magnitudes within the uncertainties. The relative contributions of the 6 Li SBU and 
DBU to the incomplete fusion/transfer process has been discussed in great details 
in Refs. Q3HQl|4r| by considering the corresponding lifetimes obtained by using 
a semi-classical approach fully described in a previous publication |[39ll . We con- 
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Fig. 5 Experimental 
(39] ESI Ell and theoretical 
CDCC (30) angular distribu- 
tions for the SBU and DBU 
projectile breakup processes 
(see text for details) obtained 
at Ei ab = 25.5 MeV and 
29.6 MeV for 6 Li+ 59 Co. The 
chosen experimental contin- 
uum excitation energy ranges 
are given. (Courtesy of F.A, 
Souza) 
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elude that the flux diverted from complete fusion to incomplete fusion would arise 
essentially from DBU processes via high-lying continuum (non-resonant) states of 
6 Li; this is due to the fact that both the SBU mechanism and the low-lying DBU 
processes from low-lying resonant 6 Li states occur at large internuclear distances 
ll39l l40l |4H . Work is in progress to study incomplete fusion for 6 Li+ 59 Co within a 
newly developed 3-dimensional classical trajectory model BTl . 



3.3 Coupled-channel calculations for reactions induced by halo 
nuclei 

As far as exotic halo projectiles are concerned we have initiated a systematic study 
of 8 B and 7 Be induced reactions data 1521 with an improved CDCC method 11481 . 
Fig. 6 displays the analysis of the elastic scattering for the 7 Be+ 58 Ni system ll52l . 
The curves correspond to CDCC calculations with (solid lines) and without (dashed 
lines) 7 Be — > a + 3 He breakup couplings to the continuum. The 6 Li and 7 Be cal- 
culations were similar, but with a finer continuum binning for 7 Be. As compared 
to 7 Be+ 58 Ni (similar to 67 Li+ 58 ' 64 Ni) the CDCC analysis of 8 B+ 58 Ni reaction ED 
while exhibiting a large breakup cross section (consistent with the systematics) is 
rather surprizing as regards the consequent weak coupling effect found to be partic- 
ularly small on the near-barrier elastic scattering. 

Recently, the scattering process of 17 F from 58 Ni target was investigated l43l 
slightly above the Coulomb barrier and total reaction cross sections were extracted 
from the Optical-Model analysis. The small enhancement as compared to the ref- 
erence (tightly bound) system 16 0+ 58 Ni is here related to the low binding energy 




Fig. 6 Ratios of the elastic scattering cross-sections to the Rutherford cross sections as a function 
of cm. angle for the 7 Be+ 58 Ni system 1521 for incident 7 Be energies of (a) 15.09 MeV, (b) 17.13 
MeV, (c) 18.53 MeV (d) 19.93 MeV and (e) 21.43 MeV. The solid and dashhed curves denote full 
and no coupling to the continuum. (This figure has been adapted from the work of Ref. [17)) 



of the F valence proton. This moderate effect is mainly triggered from a transfer 
effect, as observed for the 2«-halo 6 He HHED and the ln-halo n Be |33 in con- 
trast to the 1 p-halo 8 B+ 58 Ni reaction where strong enhancements are trigerred from 
a breakup process B31 . 
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4 Summary, conclusions and outlook 

We have investigated the fusion process (excitation functions and extracted barrier 
distributions l23l ) at near- and sub-barrier energies for the two neighbouring reac- 
tions 32 S + 90 Zr and 32 S + 96 Zr. For this purpose a new computer code named NTFus 
ll29ll has been developped by taking the coupling of the multi-neutron transfer chan- 
nels into account by using the semiclassical model of Zagrebaev flQl- 

The effect of neutron couplings provides a fair agreement with the present data 
of sub-barrier fusion for 32 S + 96 Zr. This was initially expected from the positive Q- 
values of the neutron transfers as well as from the failure of previous CC calculation 
of quasi-elastic barrier distributions without coupling of the neutron transfers ifPTl . 
With the agreement obtained by fitting the present experimental fusion excitation 
function and the CC calculation at sub-barrier energies, we conclude that the effect 
of the neutron transfers produces a rather significant enhancement of the sub-barrier 
fusion cross sections of 32 S + 96 Zr as compared to 32 S + 90 Zr. At this point we did 
not try to reproduce the details of the fine structures observed in the fusion barrier 
distributions. We believe that to achieve this final goal it will first be necessary 
to measure the neutron transfer cross sections to provide more information on the 
coupling strength of neutron transfer because its connection with fusion is not yet 
fully understood [35 1 . 

In the second part of this lecture, we have studied the breakup coupling on elas- 
tic scattering and fusion by using the CDCC approach with a particular emphasis 
on a very detailed analysis of the 6 Li+ 59 Co reaction. The CDCC formalism, with 
continuum-continuum couplings taken into account, is probably one of the most re- 
liable methods available nowadays to study reactions induced by exotic halo nuclei, 
although many of them have added complications like core excitation and three- 
body structure. The respective effects of transfer/breakup are finally outlined for 
reactions induced by l/?-halo, l«-halo and 2«-halo nuclei. 

The complexity of such reactions, where many processes compete on an equal 
footing, necessitates kinematically and spectroscopically complete measurements 
ll53l . i.e. ones in which all processes from elastic scattering to fusion are measured 
simultaneously, providing a technical challenge in the design of broad range detec- 
tion systems. A full understanding of the reaction dynamics involving couplings to 
the breakup and nucleon-transfer channels will need high-intensity RIB and precise 
measurements of elastic scattering, fusion and yields leading to the breakup itself. A 
new experimental program with SPIRAL beams and medium-mass targets is getting 
underway at GANIL. 

Acknowledgements I would like to thank A. Diaz-Torres, N. Keeley, F.A. Souza and A. Richard 
for very fruitfull discussions on many theoretical aspects of this lecture. 
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